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ABSTRACT 

The merger of a super-massive binary black hole (SBBH) is one of the most extreme events in 
' the universe with a huge amount of energy released by gravitational radiation. Although the 

I characteristic gravitational wave (GW) frequency around the merger event is far higher than 

^ ' the nHz regime optimal for pulsar timing arrays (PTAs), nonlinear GW memory might be a 

Q^l critical smoking gun of the merger event detectable with PTAs. In this paper, basic aspects 

of this interesting observation are discussed for SBBHs, and the detection numbers of their 
memory and inspiral GWs are estimated for ongoing and planned PTAs. We find that the 
expected detection number would be smaller than unity for the two-types of signals even with 
^ I the Square Kilometer Array. We also provide various scaling relations that would be useful to 

study detection probabilities of GWs from individual SBBHs with PTAs. 

_ Key words: gravitational waves — ^pulsars: general 

1 INTRODUCTION 

] A pulsar is an excellent clock in the universe, and provides us with a powerful method to directly detect gravitational waves (GWs) (Sazhin 

■ 1978; Detweiler 1979). By using multiple pulsars and suppressing noises due to independent timing fluctuations of individual pulsars, we can 
I further improve sensitivity to GWs, which generate common signals to observed pulsars (Hellings & Downs 1983). This statistical method 
' is known as a pulsar timing array (PTA), and is now considered as a promising approach to probe GWs around the nHz regime. The most 

^ ' plausible target of ongoing PTAs is the stochastic GW background made by super-massive binary black holes (SBBHs). For detecting such a 
background, a long term operation of a PTA is crucially advantageous. Since the sensitivities of ongoing PTAs have been rapidly improved, 
r> it is likely that a PTA project would succeed in detecting the background before long (see e.g. Jenet et al. 2005; Hobbs et al. 2009). 

■ Probing GWs from individual SBBHs with a PTA is also quite interesting (Lommen & Backer 2001; Sesana et al. 2009), especially in 
relation to traditional observations of electro-magnetic wave (EMW) emissions. Actually, with a constraint on the GW amplitude, a pulsar 
timing observation (PSR B 1855+09, Jenet et al. 2004) recently ruled out the previously proposed parameters of a SBBH system postulated 
from a radio observation of the galaxy 3C66B at z — 0.02 (Sudou et al. 2003). It is expected that interplays between EMW observations and 
PTAs will be biased to objects at relatively low redshift, due to observational accessibilities (see e.g. Lommen & Backer 2001). 

On the evolutionary path of a SBBH, the most violent and fascinating phase would be its final merger where a huge amount of energy 
is released by gravitational radiation. Impacts of such an extreme event would not be limited only to the GW community. For example, it has 
been actively discussed that a transient EMW signature might be associated with a merger of a SBBH (e.g. through interaction between a 
SBBH and its circumbinary disk, see Haiman et al. 2009 and references therein). In the future, there might appear a potential interpretation 
that a peculiar time-dependent phenomena observed with EMWs is related to merger of a SBBH. 

It would be very exciting to probe a merger event by observing the intense GWs with a PTA. However, the order of the characteristic 
frequency of the merger GWs is given by l/(27rAf) ~ 3 x 10"^(Af/10^A/o)"^Hz (M: the total SBBH mass), which is far higher than 
the optimal frequency regime of a PTA around ~nHz. Indeed, depending on the masses of SBBHs, GWs around the merger phase are one of 
the primary targets of the proposed Laser Interferometer Space Antenna (LISA), which is designed to have sensitivity around O.lmHz-lHz 
(Bender et al. 1997). However, LISA will not be launched before 2019. 

Fortunately, GW signals associated with the violent merger are not completely localized at the characteristic frequency 0{AI^^). 
Because of the intense and anisotropic GW emission around the merger phase, the so-called GW memory is simultaneously generated. 
At frequencies much lower than M^^, its waveform can be regarded as a (burst-like) step function profile and its non-dimensional Fourier 
amplitude he becomes independent of the frequency. Thus a GW memory might be a critical smoking gun of a SBBH merger event detectable 
with a PTA. If detected, it would surely have a broad impact on astronomy. 
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The principle aim of this paper is to discuss basic aspects on observing memory GWs from SBBHs using PTAs. For comparison, we 
also analyze inspiral GWs. The rest of this paper is organized as follows. In section 2, we provide expressions for the amplitudes of the 
inspiral and memory GWs. Then essential properties of PTA noises are mentioned in section 3. In section 4, we study the expected detection 
numbers for the memory and inspiral GWs with PTAs. In our simple formulation, we provide various scaling relations that would be useful 
to discuss PTA observation of GWs from individual SBBHs. Throughout this paper we use the geometrical units with G — c — 1. 



2 AMPLITUDES OF INSPIRAL AND MEMORY GWS 

In this section we summarize basic expressions for the amplitudes of the inspiral and memory GWs. For a SBBH (two mass; mi and m2) in 
a circular orbit, the amplitudes of the two polarization modes of inspiral GWs at a frequency / are given by the quadrupole formula (see e.g. 
Favata 2009a); 

{h+,h^), = ■' (l + cos^7,2cosJ). (1) 

Here we defined the total mass M = mi + m2, the reduced mass ratio rj = mim2/(mi + m2)^, the distance D, and the inclination angle 
I. We have ry ^ 1/4 with the equality only for mi = m2. By evaluating the energy loss through gravitational radiation, the time before 
coalescence is given by 

^^'^ = 256^^^ ^ (^) (tS^W) (i^) ^ 

Around the optimal frequency for a PTA, / ~ InHz, the time Tqw is much longer than a realistic observational period Tots = O(10)yr, 
and thus the inspiral GWs can be regarded as almost periodic signals. 

A memory GW is expressed by a net gap of the transverse-traceless components of metric, and it is generated by anisotropic energy 
emission from a source. In the present case of a SBBH, the relevant energy emission is due to gravitational radiation, and the energy carried 
by gravitational radiation is a nonlinear function of the GW amplitude. Therefore, the memory GWs of interest are often called the nonlinear 
GW memory (Christodoulou 1991; Wiseman & Will 1991; Blanchet & Damour 1992; Thome 1992; Kennefick 1994; Favata 2009a; 2009b). 
Since most of GW energy from a SBBH is emitted around the final merger phase with a time duration 0{M), the time profile of a memory 
wave can be regarded as a step function with a time resolution longer than M. The gaps of the two polarization modes for the memory GW 
are modeled by 

(/i+,/ix)™ = sin' 7(17 + cos'/, 0) (3) 

with an 0(1) parameter g determined by the history of GW emission, mostly around the merger epoch (Favata 2009b). In what follows, we 
use a re-parameterization gi2 = g/12 with gi2 ~ 1 from a recent study by Favata (2009b). 

Next we evaluate the non-dimensional characteristic GW amplitudes defined by ftc = ,/ { \ h+(J)\'^ + \hx (,f)|'}^ ■ Here we denoted 

the Fourier transformation a{f) = J^^^ a{t) exp{2inft)dt for a function a{t) in the time domain. We also introduced the notation {■ ■ ■}j 
for the angular averages with respect to the inclination /. For the inspiral signal, we have 

where the last factor (fTobsY^^ represents the effective signal amplification due to the multiple rotational cycles. In the same manner, the 
corresponding amplitude for the nonlinear memory is given by 



, gi27?M /1543 „„^,„-i6 M \ ^IGpc^ 

Thus, for a given distance D, and an observational frequency / ~ InHz, the amplitude hem for the memory GW is larger than the inspiral 
one hci in the mass range M < 10^° Af©. But this is not true at frequencies with Tqw < Tabs, as in the case of a stellar-mass black hole 
binary observed by ground based detectors. By plugging Toba ~ Tqw in eq.l?}, we obtain hd oc f^^^^ , and also have hd > hem even at 
the frequency of the last stable orbit (Kennefick 1994; Favata 2009b). 



3 NOISES OF PTAS 

In this section, we discuss the measurement noise and the background confusion noise for GW observation with a PTA. The magnitude of 
the former, Ud, relative to the non-dimensional amplitudes, he, has the simple frequency dependence cx (Rajagopal & Romani 1995). 
The power index 3/2 can be explained as follows. When we observe a pulsar, the arrival time r(t) of its pulses is modulated by a passing 
GW signal h(t) at the Earth as h{t) oc T{t). Therefore, the measurement noise nh{t) for GW observation is related to the timing noise Td{t) 
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Figure 1. Left panel: A log-log plot for the relations between the measurement noise oc f^^^ and the GW background noise cx f~^^^ for detecting 
a GW with a non-dimensional amplitude he - The frequency fx is defined as the intersection of two noises. The lowest observable frequency f^bs = '^obs 
determined by the observational time T^ta- Right panel: the total PTA noise (solid lines) before (case A) and after (case B) the detection of stochastic GW 
background. At the optimal frequency fo, the total noise level ric becomes minimum. We have fo = fobs for the case A and fo = fx for the case B. 



as nh{t) oc Td{t). Then, in the frequency domain, we have nd{f) = f oc (l'f'd(/)P)^^^. Since the timing noise rd(t) can 

be regarded as a white noise with {rdiffTdit')) cc 5{t — t'), we obtain (jfd{f)\^^ oc /^^, and finally get nd{f) oc 

For a given PTA project, the overall shape of the noise spectrum is characterized by the two parameters Tabs and ?ido- The observational 
time Tabs determines the minimum accessible frequency fobs = rj^^ and the parameter Udo fixes the noise amplitude at a pivot frequency 

fp by 

nd{f)=ndo( 4-] . (6) 



In this paper we take fp = 10~*Hz, and use the numerical values {Tobs,ndo) ~ (4.9yr, 5.1 x 10"^'') for the Parks Pulsar Timing Array 
(PPTA) and (7.7yr, 1.2 x 10^^'') for the Square Kilometer Array (SKA), extracted from figure 1 in Hobbs et al. (2009) (see also Demorest 
et al. 2009 for other projects including nanoGrav). 

Meanwhile the GW background noise nb{f) by SBBHs has a profile nb oc /^^'''^ The spectral index —2/3 can be understood with 
the definition flew oc f^nl for the normalized energy density of a GW background per logarithmic frequency interval (see e.g. Phinney 
2001). In the frequency interval, the GW energy emitted by a binary is proportional to /^''^, as derived with the Kepler's law. When we sum 
up GWs from multiple binaries, this power-law profile is unchanged, and we have 0,gw f^^^. Then we obtain Ub oc f~^^^. 

We represent the spectrum of the GW background noise by 

Ubif) = Ubpr'/^ (J-^ (7) 

with a fiducial amplitude UbF ~ 2.0 x 10"^^ and a non-dimensional scaling parameter r. Many theoretical models of structure formation 
predict that the GW background nb{f) is mainly made by the massive end of SBBHs with Al ~ 0(10'' Af©) and the parameter r would be 
in the range 0.05 < r < 20 (Jaffe & Backer 2003; Wyithe & Loeb 2003; Enoki et al. 2004; Sesana et al. 2008). Assuming the characteristic 
mass parameters M — 10^'^A^q and ri — 0.25, the comoving merger rate R is estimated as[f| 

R = Rpr = 3.7 X lO^VCpc^Vr"^ (8) 

(from an expression in Phinney 2001), corresponding to the total merger rate on our past light-cone 7^T ~ 2.2 X lO^^r yr~^ Here we 
neglected redshift dependence of the comoving rate R. 

From a continuity equation in the Fourier space, the frequency distribution of the SBBHs per comoving volume is given by ^j^y = 



^ For a given PTA project, the timing noise level / ^I'fdC/)!^^ (essentially corresponding to Strmsif) in Sesana et al. (2009)) is proportional to /^/^. 
In this paper, we examine the expected detection rates of GW signals for individual PTAs, not for a given timing noise level at each frequency (see also Fig. 2 
in Sesana et al. (2009)). 

^ While we separately use two parameters M and 'q, the chirp mass rj^^^M is more convenient for the present argument. Note also that the reduced mass 
rjM becomes maximum at r; = 1/4 for a given chirp mass. 
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R^f = ^RTgw- This expression will be used in the next section. In the same manner, the total number of the SBBHs around a frequency 
/ within the observational bandwidth fobs is given by STZtTqw fobs/iSf), and is much larger than unity at / < 10~*Hz for a plausible 
value of r. This means that an inspiral GW signal from a binary at a typical cosmological distance would be buried under the smooth GW 
background. At / > 10~*Hz, discreteness effects of SBBHs might show up for the background (Sesana et al. 2008). But, for the fiducial 
value n(,F, we have Ud > ni, for PPTA or SKA in the high frequency regime, and, in this paper, we simply neglect the discreetness effect of 
the background. 

Now we discuss the total noise spectrum nc{f) = iniLK[nd{f),nb{f)] made by the measurement noise nd{f) and the confusion 
noise nb{f) (see figure 1). The characteristic frequency fx for their intersection is solved as fx = fpnl^p^n'^Q''^^r^^^^ and nc{fx) = 
''hi^^'^hio^^^^^'^^ ■ Po'^ observational time Tabs shorter than fx^, the total noise nc{f) is determined mainly by the measurement noise 
as nc{f) — nd{f) (case A in figure 1). But, if the time Tabs is longer than and the detection of the GW background is within reach, 
the effective noise nc{f) is a piecewise power-law function (case B in figure 1). For both cases, the optimal frequency associated with the 
minimum value of the noise level nc{f) is given by fo = msLx[fx, fobs]- 



4 EXPECTED NUMBERS OF DETECTIONS 

In this section, we estimate how many detections we can expect for inspiral and memory GWs of nearby SBBHs with PTAs. The signal-to- 
noise ratio for the detection is given by SN = hc/ric from which we can inversely obtain the observable comoving distances D as well as 
the observational volumes AivD'^ /3 for a given threshold SN. Here, to deal with these geometrical quantities, cosmological effects can be 
safely neglected, since the detectable binaries would be at a relatively low redshift. 

Then, using the relevant comoving number densities (dn/dln / for the inspiral signals and RTobs for the memory signals), we obtain 
the expected numbers of detectable events in a logarithmic frequency interval as 

dN, ^ 2M^°'''Rf/''Tll^if dN„, ^ 15i3'^^ gfiM'RTobsV^ 
din/ 5i/237r2/3nc(/)3S'iV^' din/ 967680 70^/'^-K^nc{f)^SN^ 

for the inspiral (i) and the memory (m) GWs. Since we have Ud oc f '"^^^ and Ub oc f~^^^, the total detection rates A'^; and Nm are dominated 
by signals around the optimal frequency /o and we simply put A^i — {dNi / d\n f) ,andA^m ~ {dNm/d\n f) j> for our order-of-magnitude 
estimation. 

So far we have only considered SBBHs at their massive end M ~ 10^'^ A/q . But, here, we briefly comment on the mass dependence of 
the detectable binaries. For a given merger rate dR/d In M per logarithmic mass interval, the numbers of detectable binaries are proportional 
to M°' dR/ d\n M with a = 10/3 for inspiral and a = 3 for the memory signals (see eqs.l|9]l). But theoretical studies (see e.g. Sesana et al. 
2009) predict that the mass distribution dR/d In M of the merger rate would be apparently less steep than cx in the mass range lower 
than M ~ 10® M© (e.g. at M = 10* M©). Thus the detectable GW signals would be mainly made by the massive end with M > 10® M© 
and our basic prescription so far would be justified. In the following demonstration, we take M = 10^'^Mq and ri — 0.25. 

Now let us discuss the prospects for detecting inspiral and memory GWs with the two representative projects; the ongoing PPTA and 
the planned SKA. For the fiducial background at r = 1, we have the optimal point (/o,Mc(/o)) ~ (6.5 x 10~®Hz, 2.7 x 10~^^) with 
PPTA and (1.3 x 10~**Hz, 1.7 x lO"^'^) with SKA. In fi gure 2, we plot the expected event numbers as functions of the scaling parameter r 
defined relative to the fiducial background model (see eqs.l|7]l and l[8}). In this figure we set the threshold at SN = 1, since its dependence 
is straightforward. The numbers Ni and Nm change their power-law indexes at the transition point rr oc T'^^l^^^rido where the coincidence 
fx — fobs occurs for a given project. For PPTA the transition point is at tt ~ 1, since we have fx = fobs for the fiducial background level 
(consistent with figure 1 in Hobbs et al. 2009). But the point tt becomes smaller for SKA whose observational period is longer with a smaller 
measurement noise. At r < rr, we have fo = fobs and the dependencies on the rate r are very simple Ni oc rT^^^'^ and Nm oc rT^^^^, as 

o / 3/2 1/26 

the background noise is not important in this regime (case A in figure 1). At r > tt, we have Ni x r ' T^b^ and Nm oc r ' Tobs, and 
these two numbers depend very weakly on the parameter r. Interestingly, the number Nm is now a decreasing function of r. 

For each PTA project, the expected number Nm for the memory signals is 3-4 orders of magnitude smaller than the number Ni for the 
inspiral signals. As we discussed before, the comoving number densities for the merger and the inspiral signals per logarithmic frequency 
interval are given by RTobs and 8RTgw/S respectively. In the frequency regime relevant for PTAs, the former is much smaller than the 
latter, since we have Tqw 3> Tobs. Therefore, even though the detectable distance is larger for the merger signals as indicated by eqs.(4) and 
(5), their detection rate becomes smaller than that of the inspiral signals. 

Figure 2 also shows that SKA would have at least ~ 10 times larger events Ni than PPTA, but is still unlikely to detect an inspiral event 
during its operation period. These qualitative predictions for circular SBBHs would be fairly robust, and would not be changed with a more 
detailed analysis (e.g. using an elaborate mass function of SBBHs around M ~ 10® M©). 

The author is grateful to T. Tanaka for carefully reading the manuscript. He also thanks the referee for valuable comments to improve 
the draft. This work was supported by the Grants-in-Aid for Scientific Research 20740151 from the Ministry of Education, Culture, Sports, 
Science and Technology (MEXT) of Japan. 

After submission of this paper, there appeared related studies on PTAs by Pshirkov et al. (2009), and van Haasteren and Levin (2009). 
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Figure 2. The expected numbers of events Ni (solid curves) and Nm (dashed curves) with SN ^ 1 for the reference value M = 10^' ^Mq. The horizontal 
axis represents the strength r of the GW background relative to the fiducial value ni,p. The thick curves are for PPTA with the transition point at = 1 and 
the thin curves are for SKA with = 7.5 X 10^'^ We have the scaling relations Ni, Nm oc r at r < ry, and Ni oc r^/^^ and Nm oc r"^/^® at r > ry. 

They predicted higher rates for detection of the memory GWs. The difference is mainly due to the treatment of the background noise. Without 
the background noise, the detection rates are obtained by extrapolating the lines air < rr in Fig. 2 to the regime r > tt- 
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